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Abstract. Let G be an edge-weighted hypergraph on n vertices, m edges of size 
< s, where the edges have real weights in an interval [1, W]. We show that if 
we can approximate a maximum weight matching in G within factor a in time 
Tin, m, W) then we can find a matching of weight at least [a — e) times the 
maximum weight of a matching in G in time (e"^)*^'^' x 

max logn maxmi+. ..m„=m I]? T(min{n, srrij}, TTij, (e"^)*^*' '). 

In particular, if we combine our result with the recent (1 — e)-approximation algo- 
rithm for maximum weight matching in graphs due to Duan and Pettie whose time 
complexity has a poly-logarithmic dependence on W then we obtain a (1 — e)- 
approximation algorithm for maximum weight matching in graphs running in 
time (e-i)°(^)(m + n). 



1 Introduction 

A hypergraph G consists of a set V of vertices and a set of subsets of V called 
edges of G. In particular, if all the edges are of cardinality two then G is a graph. 
A matching of G is a set of pairwise non-incident edges of G. If real weights are 
assigned to the edges of G then a maximum weight matching of G is a matching 
of G whose total weight achieves the maximum. 

The problem of finding a maximum weight matching in a hypergraph is a 
fundamental generalization of that of finding maximum cardinality matching in 
a graph. The latter is one of the basic difficult combinatorial problems that still 
admit polynomial-time solutions. For hypergraphs the decision version of the 
maximum weight matching problem is NP-hard even if the edges are of size 
0(1) since it is a generalization of the problem of maximum weight indepen- 
dent set for bounded degree graphs [15]. On the other hand, polynomial-time 
algorithms yielding {d—\ + -approximation of maximum weight matching 
in hypergraphs with edges of size d are known [3]. 

The fastest known algorithms for maximum weight matching in graphs have 
substantially super-quadratic time complexity in terms of the number n of ver- 
tices of the input graph G [11, 12, 20]. For these reasons, there is a lot of interest 



in designing faster approximation algorithms for maximum weight matching 
[4-6, 14, 18, 19]. 

Recently, even fast approximation schemes for maximum weight matching 
in graphs have been presented ^ . The fastest known in the literature is due to 
Duan and Pettie [7]. It yields a (1 — e) -approximation in time 0(me~^ log^ n) 
for a connected graph on n vertices and m edges with real edge weights. The 
approximation scheme from [7] is a composition of a (1 — e)-approxiniate re- 
duction of the problem in general edge weighted graphs to that in graphs with 
small edge weights and an efficient (1 — e)-approximate algorithm for graphs 
with small edge weights. 

1.1 Our contributions 

Let G be an edge-weighted hypergraph on n vertices, m edges of size < s, 
where the edges have size real weights in an interval [1, W]. We show that if 
we can approximate a maximum weight matching in G within factor a in time 
T{n, m, W) then we can find a matching of weight at least a — e times the 
maximum weight of a matching in G in time (e~^)*^(^) x 
max log a max^^+,,,rn„=mY.iT{'^iT^{n,smj},mj,{e-^)'^'^''~ )). 

This reduction of maximum weight matching in hypergraphs with arbitraUly 
large edge weights to that in hypergraphs with small edge weights is incompa- 
rable to the aforementioned similar reduction for graphs from [7]. In particular, 
if we combine our reduction with the aforementioned (1 — e)-approximation 
algorithm for maximum weight matching in graphs from [7] whose time com- 
plexity has a poly-logarithmic dependence on W then we obtain a (1 — e)- 
approximation algorithm for maximum weight matching in graphs running in 
time {e~^)^^^\m + n). In comparison with the approximation scheme from 
[7], our approximation scheme is more truly linear in m + n, as free from the 
poly-logarithmic in n factor at the cost of larger polynomial dependence on . 

As another corollary from our approximate edge-weight reduction for hy- 
pergraphs, we obtain also some results on approximating maximum weight in- 
dependent set in graphs of bounded degree. 

1.2 Other related results 

As the problem of finding maximum weight matching in graphs is a classical 
problem in combinatorial optimization there is an extensive literature on it. It 

' In a preliminary version of this paper presented at SOFSEM Student Forum held in January 

2010 (no proceedings), an 0(n" log n) -time approximation scheme for maximum weight 
matching in bipartite graphs has been presented. 



includes such milestones as an early algorithm of Kuhn [17] just in the bipartite 
case and an algorithm of Edmond and Karp [8] running in time O(nm^), where 
n is the number of vertices and m is the number of edges in the input graph. 
Hungarian algorithm [17] can be implemented in time 0{mn + logn) with 
the help of Fibonacci heaps [9] and this upper bound can be extended to include 
general graphs [10]. 

Assuming integer edge weights in [— W, W] and RAM model with 
log(max{A^, n})-bit words, Gabow and Tarjan established 0{^/nm\og{nW)) 
and 0(V n log nm log(nW^)) time-bounds for maximum weight matching re- 
spectively in bipartite and general graphs [11, 12]. 

More recently, Sankowski designed an O(n'^M^) -time algorithm for the 
weighted matching problem in bipartite graphs with integer weights, where oj 
stands for the exponent of fast matrix multiplication known to not exceed 2.376 
[20]. His result asymptotically improved an earher upper- time bound for max- 
imum weight matching in bipartite graphs with integer weights of the form 
O(VnmVF) due to Kao [16]. 

There is also an extensive literature on fast approximation algorithms for 
maximum weight matching in graphs [4-6, 14, 18, 19]. Typically they yield an 
approximation within a constant factor between ^ and almost |, running in 
time of order m log*^^^^ n. Already the straightforward greedy approach yields 
^-approximation in time 0(m log n). 

The maximum weight matching problem in hypergraphs is known also as a 
set packing problem in combinatorial optimization [15]. By duality it is equiv- 
alent to maximum weight independent set and hence extremely hard to ap- 
proximate in polynomial time [13]. The most studied case of maximum weight 
matching in hypergraphs is that for d-uniform hypergraphs where each edge is 
of size d. Then a polynomial-time (d — 1 -|- -approximation is possible [3]. 
By duality, one obtains also a polynomial-time (d — 1 + l/d) -approximation of 
maximum weight independent set in graphs of degree d (cf. [15]). 

2 Simple edge weight transformations 

In this section, we describe two simple transformations of the edge weights 
in the input hypergraph G such that an a-approximation of maximum weight 
matching in the resulting hypergraph yields an [a — e)-approximation of max- 
imum weight matching of G. We assume w.l.o.g. throughout the paper that G 
has n vertices, m edges, and real edge weights not less than 1. The largest edge 
weight in G is denoted by W. 



Lemma 1. Suppose that there is an a-approximation algorithm for maximum 
weight matching in G running in time T{n, m, W). Then, there is an 0{n + m)- 
time transformation ofG into an isomorphic hypergraph G* with edge weights 
in the interval [1, -] such that the aforementioned algorithm run on G* yields an 
{a — e) -approximation of maximum weight matching in G in time T(n, m, j). 

Proof We may assume w.l.o.g that W > Note that the total weight of max- 
imum weight matching in G is at least W. Hence, if we transform G to a hyper- 
graph G' by raising the weight of all edges in G of weight smaller than ^ to 

TXT- 

^ then the following holds: 

1. the maximum weight of a matching in G' is not less than that in G; 

2. any matching in G' induces a matching in G whose weight is smaller by at 
most eW. 

To find an a-approximation of maximum weight matching in G', we can 
simply rescale the edge weights in G' by multiplying them by ^ . Let G* de- 
note the resulting graph. Now it is sufficient to run the asumed algorithm on G* 
to obtain an {a — f)-approximation of maximum weight matching in G. Note 
that the apphcation of the algorithm will take time T(n, m, □ 

Lemma 2. Suppose that there is an {a — e) -approximation algorithm for maxi- 
mum weight matching in G running in time T'{n' , m' , W' , e). By rounding down 
each edge weight to the nearest power ofl + e and then running the (a — e)- 
approximation algorithm on the resulting graph, we obtain an {a — 0(e))- 
approximation of maximum weight matching in G in time T'{n, m, W, e) -|- 
0{n + m). 

Proof Let e be any edge in G. Denote its weight in G by w{e) and its weight 

in the resulting graph by w'{e). We have w'{e){l + e) > w{e). Consequently, 
we obtain w'{e) > w{e) — ew'(e) > (1 — e)w(e). It follows that a maximum 
weight matching in the resulting graph has weight at least 1 — e times the weight 
of a maximum weight matching in G. Thus, if we run the asumed (a — e)- 
approximation algorithm on the resulting graph then the produced matching 
with edge weights restored back to their original values will yield an (a — 2e)- 
approximation. □ 

3 A transformation into an (a — e) -approximation algorithm 

A subhypergraph of a hypergraph H is any hypergraph that can be obtained 
from H by deleting some vertices and some edges. A class C of hypergraphs 



such that any subhypergraph of a hypergraph in C also belongs to C is called 
hereditary. 

In this section, we present a transformation of a hypothetic a-approximation 
algorithm for maximum weight matching in a hereditary family of hypergraphs 
with edges of size 0(1) into a (a — e)-approximation algorithm. The run- 
ning time of the (a — e)-approximation algorithm is close to that of the a- 
approximation algorithm in case the largest edge weight is e^*^*^*^ •* . 

Theorem 1. Suppose that there is an algorithm for a maximum weight match- 
ing in any hypergraph having edges of size < s and belonging to the same 
hereditary class as G running in time T{n' , m' , W) = Q{n' + m'), where n', 
ml are respectively the number of vertices and edges, and [1, W'] is the interval 
to which all edge weights belong. There is an {a — e) -approximation algorithm 
for a maximum weight matching in G running in time (e^^)*^^^^ x 
max log^ , maxn,,+...^^=n, X)i T{mm{n, smj},mj, (e~^)'^("~')). 

Proof. We may assume w.l.o.g that W = 0{n/e) and any edge weight is a 
nonnegative integer power of 1 + e by Lemmata 1, 2. Order the values of the 
edge weights in G in the increasing order. Set k = 0(e^^) and I = [log]^.,.^ -] . 
By the form of the edge weights and the setting of /, the following holds. 

Remark 1: For any two different edge weights wi and W2, if the number of 
wi is greater than that of W2 by at least I in the aforementioned ordering then 

^Wi > W2. 

In order to specify our (q — e)-approximation algorithm, we partition the 
ordered edge weights into consecutive closed basic intervals, each but perhaps 
for the last, containing exactly / consecutive edge weights, see Fig. 1. 
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Fig. 1. Partitioning of edge weights (/ = 3) 



Next, we group A;-tuples of consecutive basic intervals into large intervals 
composed of — 1 consecutive basic intervals followed by a single basic interval 
called a gap. This partition corresponds to the situation when the so called shift 
parameter x is set to 0. For x G {1, .., k — 1}, the partition into alternating large 



intervals and gaps is shifted by x basic intervals from the right, so the first large 
interval from the right is composed solely ofk — 1 — x basic intervals, see Fig. 
2. The maximal subgraph of G containing solely edges in the large intervals in 
the partition is denoted by Gx ■ 
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Fig. 2. An example of shift: 1=3, k=3 
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For our (a — e) -approximation algorithm for a maximum weight matching 
in G see Fig. 3. We shall assume the definitions of the subgraphs G^, Gxj, 
given the algorithm. 

Since the union of the gaps over all shifts covers all weights there must a 
shift where the gaps cover at most i of the weight of optimal matching of G. 
Hence, there must be a shift x such that the weight of optimal matching in Gx is 
at least (1 — l//c) of the weight of optimal matching of G. Thus, it is sufficient 
to show that Mx closely approximates an a-approximate weight matching of 
Gx- 

Consider a maximum weight matching OMx of Gx and the a-approximation 
Mx.j of a maximum weight matching of Gx.j, respectively. Note that Mxj has 
total weight not smaller than a times the total weight of OMx restricted to the 
edges in Gxj- On the other hand, each edge e in Mxj can eliminate at most 
0(1) edges of OMx from all Gx,i for i > j. The total weight of the at most 
0(1) edges is only at most the e fraction of the weight of e by Remark 1. Let 
EOMx denote the set of all edges in OMx eliminated by Mx = Uj ^x,j- The 
following two inequalities follow: 



Algorithm 1 



1. for a; -f— 1 to A: — 1 do 

2. <- dl; 

3. ^ G,; 

4. fori^ ltoO(logi+,f)do 

5. begin 

6. Set Gx,j to the sub-hypergraph of G^ induced by the edges whose weights 

7. fall in the jth interval from the right; 

8. Run the a-approximation algorithm for maximum weight matching M^j of G; 

9. ^ M^lJM^j; 

10. Remove all edges incident to M^j from G^; 

11. end 

12. Return the heaviest among the matchings 

Fig. 3. The (a — e) -approximation algorithm. 



weight{Mx) + weight{EOMx) >ax weight{OMx) 
€ X weight{Mx) > weight{EOMx) 
Consequently, we obtain: 

weight{Mx) > axweight{OMx)—exweight{Mx) > {a—e)xweight{OMx) 

Thus, Mx approximates within (a — e) a maximum weight matching of 
Gx, and consequently the heaviest of the matchings Mx approximates within 
(1 — e)(l — 1/fe) a maximum weight matching of G. By setting k = we 
obtain an (1 — 0(e)) -approximation of the optimum. 

It remains to estimate the time complexity of our method. Note that the 
weight of heaviest edge in Gx,j is at most 

+ = 0(e-i)0(e-) = (e-i)0(^-^) 

times larger than that of the lightest one. Let mx^j denote the number of edges 
in Gx,j ■ Next, let nx^j denote the number of vertices in the sub-hypergraph of 
Gx,j induced by the edges of Gxj- Note that nx,j < mm{n,smxj} by our 
assumption on the size of edges in G. 

Hence, by rescaling the weights in Gxj, we can find Mxj in time 
T(m.m{n, smxj},mxj, (e^^)'^^'^ for j = l..,0{logi_^^j/lk) and x = 
0, k-1. Note that logi+, f = = ©(e'^ log f ) and similarly Ik = 



logi+, f 0(e-^) = 0(i^^e-^) = e{e-^ loge"^). It follows that for a given 

X, the largest value of j, i.e., the number of the subgraphs Gxj is 0{e^^). 

Note that m^.j < m since each edge of G belongs to at most one hyper- 
log — 

graph Gxj- Thus, the total time taken by finding all j for j = 1, 0{ej^^^) 
for a fixed x is 

log c ^ 

Recall that x ranges over 0{e~^) possible values. 

By the assumed form of the edge weights in G, we can apply a standard 
radix sort with 0(e~^ log ^) buckets to sort the edges of G by their weights in 
time 0(m + e"^ log f ). The latter is also 0{€-'^T{n, m, (e"^)*^^^"'))) by the 
assumptions on T. 

In order to efficiently construct the graphs Gx,j, the sorted edge list is kept 
in array and there are double Unks between an occurrence of an edge in the 
adjacency lists representing G and its occurrence in the sorted edge list. To de- 
termine the edges inducing G^j, we just scan a consecutive fragment of the 
sorted list from left to right. Given a list of edges of Gx,j, an adjacency rep- 
resentation of the sub-hypergraph can be constructed in time 0{n + m) = 
0(r(n,m,(e-i)0(^"')))by using the aforementioned double Unks. 

To remove an edge from G'^ , we locate it on the sorted edge list by using the 
double hnks with the adjacency hsts and then link its predecessor with its suc- 
cessor on the sorted list. We conclude that the updates of G'^ take time 0{m) = 
0(r(n,m,(e-i)0(^"'))). □ 



4 Applications 

There are at least two known exact algorithms for maximum weight matching 
in bipartite graphs with integer edge weights for which the upper time bounds 
on their running time in linear fashion depend on the maximum edge weight W 
[16, 20]. Recently, Duan and Pettie have provided substantially more efficient 
1 — e approximation algorithm for maximum weight matching in general graphs 
with integer edge weights, whose running time also depends on W in linear 
fashion [7]. Furthermore, their final approximation scheme for this problem in 
fact exhibits poly-logarithmic dependence on W. 

Fact 1 (Duan and Pettie, see the proof of Theorem 1 in [7]). An (1 — e)- 

approximation of maximum weight matching in a connected graph on m, edges 
and positive integer weights not exceeding W can be found deterministically in 
time 0(e~^m log^ W). 



We can trivially generalize the upper time bound of Fact 1 to include a non- 
necessarily connected graph by extending it by an additive factor of 0(n). 

There is one technical difficulty in combining Facts 1 with Theorem 1. 
Namely, in the theorem we assume that there is available an a-approximation al- 
gorithm for maximum weight matching for graphs belonging to the same hered- 
itary class as G with arbitrary real edge weights not less than 1 whereas the 
algorithm of Facts 1 assumes integer weights. In fact, even if the input graph 
got positive integer weights the preliminary edge weight transformations in the 
proof of Theorem 1 would result in rational edge weights. There is a simple rem- 
edy for this. We may assume w.l.o.g that e is an inverse of a positive integer and 
through all the steps of our approximation scheme round down the edge weights 
to the nearest fraction with denominator 0(e~^) and then multiply them by the 
common denominator to get integer weights. This will increase the maximum 
weight solely by 0(e~^) and will preserve close approximability. 

Hence, Fact 1 combined in this way with Theorem 1 yield our main appli- 
cation result by straightforward calculations. 

Theorem 2. There is an approximation scheme for a maximum weight match- 
ing in a graph on n vertices and m edges running in time {e~^)^^^\m + n). 

5 Extensions 

Note that Theorem 1 includes as a special case the problem of finding a max- 
imum weight independent set in a graph G of maximum degree d which is 
equivalent to the problem of finding a maximum weight matching in the dual 
hypergraph with edges corresponding to the vertices of G and vice versa. 

Several combinatorial algorithms for maximum independent set achieving 
the approximation ratio of 0{d), where d is the maximum or average degree are 
known in the literature [15]. In the appendix, we demonstrate that by using the 
method of Theorem 1 they can be simply transformed into good approximation 
algorithms for maximum weight independent set. 
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7 Appendix: Approximation algorithms for maximum weight 
independent set in bounded degree graphs 

Note that Theorem 1 includes as a special case the problem of finding a max- 
imum weight independent set in a graph G of maximum degree d which is 
equivalent to the problem of finding a maximum weight matching in the dual 
hypergraph with edges corresponding to the vertices of G and vice versa. 

Several combinatorial algorithms for maximum independent set achieving 
the approximation ratio ofO{d), where d is the maximum or average degree are 
known in the Uterature [15]. Here, we demonstrate that by using the method of 
Theorem 1 they can be simply transformed into good approximation algorithms 
for maximum weight independent set. 

Lemma 3. Suppose that there is an a{d)-approximation algorithm for maxi- 
mum independent set in a graph on n vertices and maximum ( or average degree, 
respectively) degree d running in time S(n,d), where the function S is non- 
decreasing in both arguments. There is an a{dW) -approximation algorithm for 
maximum weight independent set in a graph on n vertices, maximum ( or av- 
erage degree, respectively) degree d, positive integer weights not exceeding an 
integer W, running in time S{nW, dW). 

Proof: Let G be the input vertex weighted graph G. We form the auxiliary 
unweighted graph G* on the base of G as follows. In G*, we replace each vertex 
V in G with the number of its copies equal to the weight of v. We connect 
each copy of v by an edge with each copy of each neighbor of v. Next, we run 
the assumed algorithm for maximum unweighted independent set on G* . Note 
that any maximal independent in G* is in one-to-one correspondence with an 
independent set in G since whenever a copy of v is in the independent set then 
all other copies of v can be inserted into it without any conflicts. □ 

The drawback of Lemma 3 is that the approximation factor and/or the run- 
ning time of the resulting algorithm for the weighted case can be very large 
in case the maximum weight W is large. However, we can plug Lemma 3 in 
the method of Theorem 1 to obtain much more interesting approximation algo- 
rithms in the weighted case. 

Theorem 3. Suppose that there is an a{d)-approximation algorithm for maxi- 
mum independent set in a graph on n vertices and maximum degree d running 
in time S{n, d), where the function S is non-decreasing in both arguments and 
S{n,d) = f2(nd\ogn). There is an {a{d€~^)'~'^^ ^) — de)-approximation al- 
gorithm for maximum weight independent set in a graph on n vertices, with 
maximum degree d, positive integer vertex weights, running in time 
0(eig^5(n(6-i)0(-^),d(6-i)«(-^))). 



Proof, sketch. Recall that the problem of maximum (weighted or unweighted) 
independent set is equivalent to the problem of maximum (weighted or un- 
weighted, respectively) matching in the dual hypergraph. In the dual hyper- 
graph, the edges have size not exceeding the maximum vertex degree in the 
input graph. We run the method of Theorem 1 on the dual hypergraph using as 
the black box algorithm the result of the application of Lemma 3 to the assumed 
algorithm and its adaptation to the maximum matching problem in the dual hy- 
pergraph. □ 



